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Lattice analysis of semi-leptonic form factors* 
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We present preliminary results from simulations done on 170 32'' x 64 lattices at f3 — 6.0 using quenched Wilson 
fermions. This talk focuses on the behavior of the form-factors, extrapolation in quark masses, dependence 
on renormalization scheme, and comparison with heavy-quark effective theory (HQET). Even though we cannot 
estimate errors due to quenching and discretization, our results are consistent with experimental results for D 
decays. We present results for the Isgur-Wise function and estimate (,'{w = 1) = 0.97(6). 



1. TECHNICAL DETAILS 

We briefly mention some details of our analysis 
to extract the form-factors at = relevant to 
phenomenology. A full analysis will be presented 
elsewhere The details of the lattices are con- 
tained in the papers on the hadron spectrum 
H), and decay constants Preliminary results 
on a sub-set of lattices have been presented at 
LATTICE 94 [| and DPF 94 §. We do not have 
data to extrapolate to nib or to a = 0, thus our 
results are relevant for D decays, i.e. D — Kliy, 
D nil/, D — > K*lv, and D plv, calculated 
at a-^ = 2.33(4) GeV. 

The decaying D meson is created at rest by 
using a p = source. On each lattice we make 
two measurements by creating the D meson at 
t = 7 and 57, for a total statistical sample size 
of 340. The sink for the final state meson is at 
t = 32 in both cases. The time-slice of the weak 
operator is taken to vary between 10 < t < 30 
and 35 < i < 55 respectively, and the insertion 
is at 5 lowest values of lattice momenta. The 
quark propagators are created using a Wuppertal 
smeared source as described in 

In order to isolate the desired matrix element 
(ME) we construct a ratio of 3-point to 2-point 
correlation functions. We have a choice of using 
either smeared-smeared {SS) or smeared-local 
{SL) 2-point correlation functions. We calculate 



*Based on talks presented by Tanmoy Bhattacharya and 
Raj an Gupta. These calculations have been done on the 
CMS at LANL as part of the DOE HPCC Grand Challenge 
program, and at NCSA under a Metacenter allocation. 



Figure 1. Ratio of 3-point to 2-point SL function 
versus the time slice of the operator insertion. 
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the AIE both ways and take the average as our 
best estimate. Figure ^ shows a typical exam- 
ple of signal for the ratio of correlators in the 
SL case: the quality in the SS case is very sim- 
ilar. We have seen a steady improvement in the 
consistency between these two estimates of ME 
with statistics. With the current sample they are 
within 1(7 in all cases. 

Pole dominance hypothesis (PDH): It states 
that all form-factors, /(Q^), have the structure 



/(Q') - f{0)/{i-QVM'), 
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where AI is the mass of the nearest resonance 
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Figure 2. Three fits to test behavior of ]^ 
1. 
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with the right quantum numbers. To test PDH 
we make two kinds of fits: (i) single parameter 
"pole" fit where M is the lattice measured value 
of the resonance mass, (ii) two parameter "best" 
fit where M and /(O) are free parameters. An 
example of these fits is shown in figure |. Overall 
we find that only /q and jv are well described 
by the "pole" form. and /^(, are consistent 
with "pole" form with M < Mpoie, while for 
we find M > Mpoie- /ai and show a much 
smaller dependence than expected from pole 
dominance, however the data are too noisy to 
make a definite statement. We use results from 
"best" fit for our final estimates. 
HQET: At leading order in l/rric, HQET pre- 
dicts that the Isgur-Wise function ^ describes all 
form-factors relevant io D ^ Klv and D K*lv 
decays. Neglecting 0{as) corrections, one gets 
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where Mi and Mf are the initial and final meson 
masses and 
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Figure 3. Comparison of f+{q^) (squares) and 



(1 - (m ImA'^ ) ^foil^) (diamonds) versus w. 
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The first HQET relation, (1 - q'^/{Mi + 
-^^/)^)^^/o = /+! is satisfied by our data as exem- 
plified in figure H for CC/3 U1U3 decay. Thus, 
both /o and /+ cannot simultaneously satisfy the 
PDH. If leading order HQET holds, i.e. Mi + 
Mf ^ mc + iris Mpoie, f+ should obey a 'dipole' 
form as fo agrees with pole dominance. The data, 
as shown in figure ||, suggests that /+ lies in be- 
tween "pole" and "dipole" forms. Similar analysis 
for the vector form factors is under progress. 
Dependence on quark mass: Figure ^ shows 
an example of the variation of /+ (0) with quark 
masses. There is significant dependence on the 
mass of the quark C decays into. This is a kine- 
matic effect as shown in Section ^ Our data 
is good enough to expose slight dependence on 
TTispectator — the Small dccrcasc in slope between 
the transitions CUi — > UiUi and CUi SUi is 
consistent with HQET. 

Fixing ms'. We fix mg using Mj./Mp and 
M^/Mp. Our preferred way is ms{M^), and the 
variation of results with is shown in Table |^. 
0{ma) effects: In Ref. Q we discuss 0{ma) 
effects for heavy-light mesons. The kinetic mass 
M2 = dPE/dp^ is given by sinhM as the data 
agree with the dispersion relation sinh^i?/2 = 
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Figure 4. Extrapolation of f+{Q^ — 0) to m„. 
f+{D Klv) is extracted from points labeled 
by squares and octagons, while f+{D ttIv) is 
from data with degenerate qq points (crosses) . 




0.01 0.02 0.03 0.04 



Table 1. Estimates of form factors in 3 commonly 
used renormalization schemes defined in Vm . 





TADl 


TADtt 


TADUo 


/+ 


0.71(4) 


0.75(4) 


0.66(4) 


/o 


0.73(3) 


0.77(3) 


0.68(2) 


fv 


1.28(7) 


1.36(7) 


1.19(6) 


fAo 


0.84(3) 


0.85(3) 


0.79(3) 


fAi 


0.72(3) 


0.74(3) 


0.68(3) 


fA2 


0.49(9) 


0.50(9) 


0.46(8) 


JAs 


0.85(3) 


0.87(3) 


0.80(3) 


fv/fAi 


1.78(7) 


1.84(8) 


1.76(7) 



sin^ p/2 + sinh^ M/2. We show variation of form- 
factors with M in Table ^. 

Renormalization Constants: To relate lattice 
results to experimental data we need the renor- 
malization constants Zj^ and Zy. We use three 
Lepage-Mackenzie tadpole improved schemes de- 
scribed in Ref. Q]. Our preferred scheme is 
TADl, and the variation with the schemes is il- 
lustrated in Table |l|. 

2. RESULTS at /3 = 6.0 

Our final results for D ~* Klv and D K*lv 
are shown in tables |l| and |^ along with vari- 
ation with TOs, type of fit, heavy-light meson 



Table 2. Estimates in TADl scheme at = 0. 
The variations give estimates of systematic errors. 







Pole 


Best 






ms[MKj 


ms(M^) 


vis [Mk j 


ms{M^) 


J+ 


Ml 


0.78(2) 
0. (9(2) 


0.80(2) 

0.<3U(z ) 


0.70(4) 

0.0((5j 


0.71(4) 
0.6o(4j 


Jo 


Ml 


n 7(^0^ 

U. I U^Z J 

0.72(2) 


n 70(9) 

0.73(2) 


68(^31 


0.70(2) 


J V 


Ml 

M2 


1.27(7) 
1.22(6) 


1.28(6) 
1.23(6) 


1.27(7) 
1.22(9) 


1.28(7) 
1.23(9) 


fA„ 


Ml 

M2 


0.85(3) 
0.84(3) 


0.85(3) 
0.84(3) 


0.83(3) 
0.80(3) 


0.84(3) 
0.81(3) 


fAi 


Ml 
M2 


0.67(2) 
0.65(2) 


0.68(2) 
0.66(2) 


0.71(3) 
0.67(4) 


0.72(3) 
0.69(4) 


fA2 


Ml 
M2 


0.46(9) 
0.42(9) 


0.47(8) 
0.44(9) 


0.49(9) 
0.50(14) 


0.49(9) 
0.52(12) 


JAs 


Ml 
M2 


0.84(3) 
0.81(3) 


0.85(3) 
0.82(3) 


0.84(3) 
0.81(4) 


0.85(3) 
0.82(4) 


h 
U 


Ml 


0.93(1) 


0.94(1) 


1.02(4) 


1.02(3) 


M2 


0.95(1) 


0.95(1) 


1.03(4) 


1.02(4) 


iv 


Ml 


1.83(7) 


1.81(6) 


1.79(8) 


1.78(7) 




M2 


1.85(7) 


1.83(6) 


1.81(11) 


1.79(10) 


iA2 
fAi 


Ml 


0.70(10) 


0.70(9) 


0.69(12) 


0.68(11) 




M2 


0.65(11) 


0.67(11) 


0.74(18) 


0.75(15) 



mass, and the renormalization scheme. Our pre- 
ferred estimates are with ms{M^), "best" fit. Mi, 
and TADl scheme. We also get that fl/f^ = 
0.87(4). 



3. dr(g2) 

As explained in section ^ there are considerable 
uncertainties involved in extrapolating the form 
factors to = 0. Therefore, we also calculate 
dT{Q^) / dQ^ by linearly extrapolating in the 
form-factors at fixed 3-momentum transfer. In 
figure H we show the results for D Keve- In fig- 
ure ^ we show the longitudinal, transverse and the 
total decay widths for the process D — > K*eVf,. 

4. Isgur-Wise function 

As mentioned in section at the leading order 
in the heavy quark mass, there is one universal 
Tsgur-Wise' function ^ which controls all the form 
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Figure 5. {l/V^^)dT{q^)/dq^ versus in GeV^ 
The shape is in qualitative agreement with exper- 
imental data. 
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factors. In figure ^ we present data for ^ren(w), 

h+{w) = [Ci{w) + ^{C2{w) + C^{W))\ ireniw) 

where Ci, C2, 6*3 are HQET rcnormalization con- 
stants including the leading 0{as) corrections 
0. The data show a slight dependence on the 
spectator quark mass indicative of correc- 
tions. Analysis from the vector form factors is in 
progress. From this data we estimate the slope 
to be Ciw = 1) = 0.97(6). 
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Figure 6. {l/V^g)dr{q^)/dq^ versus in GeV 
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Figure 7. (,{w)/^{l) at various quark masses and 
momenta. The symbols label the flavor that the 
C quark decays to, and variation of ^ with w is 
a kinematic effect. Dependence on m spectator is 
shown most clearly by the clusters of 3 points 
at w K, 1.2. The largest value in each cluster 
corresponds to the lightest spectator {U^). The 
data are fit to £,{w) = ^exp(-(2p2 _ 1)h^), 
from which we estimate the slope 
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